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ABSTRACT 

Context. The Kepler mission has provided high-accurate photometric data in a long timespan for more than two hundred thousands 
stars, looking for planetary transits. Among the thousands of candidates detected, the planetary nature of around 15% has been 
established or validated by different techniques. But additional data is needed to characterize the rest of the candidates and reject 
other possible configurations. 

Aims. We started a follow-up program to validate, confirm, and characterize some of the planet candidates. In this paper we present 
the radial velocity analysis of those presenting large variations, compatible with being eclipsing binaries. We also study those showing 
large rotational velocities, which prevents us from obtaining the necessary precision to detect planetary-like objects. 

Methods. We present new radial velocity results for 13 Kepler objects of interest (KOIs) obtained with the CAFE spectrograph at the 
Calar Alto Observatory, and analyze their high-spatial resolution (lucky) images obtained with AstraLux and the Kepler light curves 
of some interesting cases. 

Results. We have found five spectroscopic and eclipsing binaries (group A). Among them, the case of KOI-3853 is of particular 
interest. This system is a new example of the so-called heartbeat stars, showing dynamic tidal distortions in the Kepler light curve. 
We have also detected duration and depth variations of the eclipse. We suggest possible scenarios to explain such effect, including the 
presence of a third substellar body possibly detected in our radial velocity analysis. We also provide upper mass limits to the transiting 
companions of other six KOIs with large rotational velocities (group B). This property prevents the radial velocity method to obtain 
the necessary precision to detect planetary-like masses. Finally, we analyze the large radial velocity variations of other two KOIs, 
incompatible with the presence of planetary-mass objects (group C). These objects are likely to be stellar binaries. However, a longer 
timespan is needed to complete their characterization. 

Key words. Techniques: radial velocities, Planets and satellites: general, (Stars:) binaries: eclipsing, close. 


1. Introduction 


The need of ground-based follow-up observations of the planet 
candidates provided by the Kepler mission dBorucki et al.ll2010l) 
was known even prior to its launch on March 2009. This need 
comes from two sides. First, high spatial resolution is required 
owing to the large pixel size (~ 4 arcsec) and the typical photo¬ 
metric apertures (6-10 arcsec) used by Kepler, which permits the 
contamination of light-curve by close sources, micking planet¬ 
like transits. Second, radial velocity (RV) observations of the iso¬ 
lated candidates are crucial to identify or reject their planetary 
nature. The Community Follow-up Observing Program (CFOFLji 
serves as a communication tool between the different interna¬ 
tional groups working on this. The false positive rate in the Ke¬ 
pler sample of pl anet candidates has been investigated by sev- 
eral authors (e.g., iMorton & Johnsonl201 lllSanterne et al .120 1 2t 


iFressin et al.l 120131: ICoughlin et al.ll2014l) . and more than 3000 
Kepler Objects of Interest (KOIs) have been rejected so far by 
different groups using different techniques^. 

In 2011, we started a ground-based follow-up pro¬ 
gram with both high-resolu t ion (lu cky) imaging with As¬ 
traLux dLillo-Box et al] 120121 l2QI4al) to detect possible close 
companions and hig h-resolution spectroscopy with CAFE0 
(lAceituno et al.l 120131) to characterize some of the Kepler can¬ 
didates via RV studies, both instruments installed at the 2.2m 
telescope in the Calar Alto Observatory. 

In this paper, we present the results of part of our high- 
resolution spectroscopic survey of isolated planet candidates 
in the Kepler sample. This survey has already confirmed the 
planetary nature of Kepler-9 lb (KO I-2133b. iLillo-Box et aD 
l2014dbh . Kepler-432b ( KOI-1299b, ICiceri et al.l 1201 5l) . and 
Kepler-447b (KOI-1800b. ILillo-Box et al.l 120151) : and has de- 


* Based on observations collected at the German-Spanish Astronom¬ 
ical Center, Calar Alto, jointly operated by the Max- Planck-Institut fur 
Astronomie (Heidelberg) and the Instituto de Astroffsica de Andalucfa 
(IAA-CSIC, Granada). 

1 https://cfop.ipac.caltech.edu/home/ 


2 According to the CFOP and the NASA Exoplanet Archive 
(http://exoplanetarchive.ipac.caltech.edu lAkeson et al.l 12013 ). among 
the 7348 identified KOIs to date, 3170 have been labelled as "False pos¬ 
itives". 

3 Calar Alto Fiber-fed Echelle spectrograph 
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tected the secondary eclipse of WASP-10b (ICruz et al.ll2015l) . 
In the present paper, we analyze 13 specific candidates in our 
sample that have shown either large RV variations or broad¬ 
ened spectra owing to the fast rotation of the host stars. Among 
the studied candidates, we have found an eclipsing binary that 
shows heartbeat-like modulations on it_s light curve. This effect 
was explained and explored by iKumar et al.l (! 1995 1) and is due 
to tidal distortions in eccentric binary systems. The first object 
discovered whose light curve is m a inl y dominated by this type 
of modulation was KOI-54 ( IWelsh et al.ll201 ll) . Since then, only 
a few more cases have been discovered (e.g., [Thom pso n et al .l 
l2012h . but more cases are needed to understand and test theo¬ 
retical models of how pulsations are induced from the dynamic 
tidal distortion. Interestingly, the RV analysis of KOI-3853 also 
shows indications of a possible third companion. 

This paper is organized as follows. In Sect. 32.31 we describe 
the observations used in this analysis. Section §j3]shows the anal¬ 
ysis and results of the RV data for the found undiluted binaries 
(Sect. 33.1b . Also, we set upper limit masses for the companions 
of few fast rotators (Sect. 33. 2b and study two yet unsolved cases 
with current data in Sect. 33.31 In Sect. ^4]we show the detection 
of peculiar (not yet reported) features in the light-curves of some 
of the targets analyzed in this paper. Finally, a brief discussion 
and conclusions of these findings is provided in Sect. |5] 

2. Observations 

2.1. Kepler photometry 

The objects analyzed in the present paper were firstly classified 
by the Kepler mission as potential planet hosts by the detection 
of planetary-like transits. We have used the long cadence data 
of quarters Q1-Q17 of the Kepler mission to analyze the light 
curves (LCs) of some of these Kepler objects of interest (here¬ 
after KOIs). The data were obtained from the MASlQ (Milkuski 
Archive for Space Telescopes) database. We used the Pre-search 
Data Conditioning Simple Aperture P hotometry flux (PDCSAP 
ISmith et all2012nStumne et al.il2012h and detrend these data by 
using simple spline fitting to every set of data without jumps 
longer than one day. We have used one-day averaged flux bins 
as nodes for the spline fitting. In specified cases (such in KOI- 
3853), we removed some parts of the light-curve in the calcula¬ 
tion of the nodes in order to not smooth or even eliminate non 
artificial modulations (such as primary transits, eclipses or other 
real variations). 

2.2. High-spatial resolution images 

The presence of undiluted sources can contaminate both the Ke¬ 
pler light curve and the spectra of the target, introducing imprints 
that could mask or mimic a planetary signal or conta minat e 
the determination of its properties (e.g. iDaemgen et alJl2009t 
lAdams et al .11201 2l lLaw et alll20l4l) . Thus, prior to the RV ob¬ 
servations, we obtained high-resolution images for the targets by 
using the lucky-imaging technique with the A straLux instrument 
at the 2.2m telescope in Calar Alto. In iLillo-Box et al.l (120121 
'2014c ) we already published the results of our high-spatial reso¬ 
lution survey of the Kepler candidates and provided details on 
the reduction and analysis of the data. More than 170 KOIs 
were observed and analyzed, including 9 out of the 13 KOIs 
studied here. For the re main ing target^ we used ancillary high- 
resolution images from lLaw et al.l(l2014l) and the UKIRT J-band 

4 https://archive.stsci.edu/kepler/data_search/search.php 


Table 1. Ancillary data and calculated BSC (blended source confidence) 
values for the KOIs studied in this paper. 


KOI High-Res lm. a Comp. Sep. Am BSC 

Ref. (< 6") arcsec. mag (band) % 


Group A - Eclipsing binaries 


340.01 

L13 

0 h 

99.99 

371.01 

L13 

0 b 

96.0 

686.01 

LB 14 

0 

99.9 

3725.01 

LB 14 

0 

99.7 

3853.01 

LB 14 

0 

99.6 


Group B - Fast rotators 


12.01 

LB 14 

0 

99.3 

131.01 

LB 12 

1 

5.54" 3.48 (isDss) 

366.01 

LB 14 

0 

99.3 

625.01 

LB 14 

0 

96.8 

972.01 

CFOP 

0 

- 

3728.01 

LB 14 

0 

99.6 


Group C - Unsolved cases 

1463.01 CFOP 

1 

5.47" 

6.9 (Kep.) 

- 

3890.01 LB 14 

0 



98.3 


Notes. Uncertainties in the angular separation of the companions and in 
their magnitude difference have been omitted for clarification purposes. 
We refer to the corresponding referen ces for more informatio n about 
t his. (n) LB 12 =lLillo-Box et alldMl) . L13 = lLaw et al.ld20li) . LB 14 
- ILillo-Box et alJl2014ct) rCFQP = no reference found in the bibliogra¬ 
phy but companion tables of the UKIRT J-band survey by D. Ciardi are 
availa ble in the CFOR The maximum angular separation explored in 
Law_et al.j d20ll is 2.5 arcsec. 

surve\0 of the Kepler sample, publicly available in the CFOP 
and provided by D. Ciardi. 

In Table [T| we summarize the main results obtained for each 
of the KOIs of the present study, including the calculated blended 
source confidence (BSC) and the number of detected visual com¬ 
panions closer than 6.0 arcsec. The BSC parameter was pre¬ 
sented in lLillo-Box et al .l (12014a i and represents a standard and 
unbiased procedure to derive the probability of a given source to 
be isolated (i.e. without any blended undetected sources). Addi¬ 
tional sources were found closer then 6.0 arcseconds in just two 
KOIs (KOI-0131 and KOI-1463). However, these visual com¬ 
panions are too faint to affect the CAFE spectrum and thus their 
contribution to the RV is negligible. Added to this, they lie com¬ 
pletely outside of the CAFE fiber (with 2.1 arcsec of diameter), 
thus not contaminating the spectra of the KOIs. 

2.3. CAFE observations and data reduction 

We have obtained multi-epoch RV data for the 1 3 KOIs shown 
in Table Q] We used the CAFE instrument dAceituno et al.l 
120131) mounted on the 2.2m telescope at Calar Alto Observa¬ 
tory (Almeria, Spain) to obtain high-resolution spectra (R = 
59000-67000) in the optical range (4000-9000 A). Table[2]sum- 
marizes the targets observed, the number of epochs acquired, and 
the time span (defined as the time between first and last spectrum 
obtained) for each of them. 

We have used the standard reduction process for the 
CAFE spectra, already described in lLillo-Box et al.l (1201 4bl) and 
ILillo-Box et al.l d2014cl) . The basic reduction is performed by us- 


5 http://keplergo.arc.nasa.gov/ToolsUKIRT.shtml 
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ing the observatory pipeline described in lAceituno et alJ (1201 310 . 
The wavelength-calibrated spectra are then analyzed with our 
own routines in order to obtain the RV. In the case of the rapid 
rotator s, we_used the genetic algorithm GAbox (presented in 
[Lillo-Box et al.l [2014cill to determine the RV at the different 
epochs. We used the theoretical equivalency between maximum- 
likelihood and cross-correlation presented in IZuckerl (120031) to 
match the observed spectrum to a synthetic template of the same 
spectral type. T his RV extraction process is fully described in 
iLillo-Box et~aO d2014bl) . 

For the slow rotators, a more precise anal ysis c an be done 


by using the cross-correlation technique dBaranne et a l. 


We used a binary weighted mask as template dPepe et al 


1996i). 


20021 ) 


and cross-correlate the observed spectra against it. The mask was 
designed by using a solar spectrum obtained from the solar at¬ 
las provided by BASS 200(0- From this spectrum, we selected 
a series of isolated, sharp, and high-contrast absorption lines. In 
total, our mask includes around 2100 lines in the CAFE spectral 
range. We used a velocity range of +30 km/s around the expected 
RV value for each case. The RV of each spectrum is measured as 
the peak of a gaussian fit to the resulting cross-correlation func¬ 
tion (CCF). This velocity is then corrected from the barycentric 
Earth radial velocity (BERV) obtained from the corresponding 
Julian date at mid-exposure. 

In Tables lAJ~l to lA.13l of Appendix lAl we provide the calcu¬ 
lated radial velocities and the observing information (signal-to- 
noise, julian date, and phase) for each epoch and KOI. 


3. Radial velocity analysis 

The extracted RV data explained in the previous section was an- 
alyzed by using the RVLIN softwai'tQ dWright & Howardll2009b 
and its additional package BOOTTRAN for param eter uncer¬ 
tainties estimation with bootstrapping dWang et al.ll2012 ) to fit 
our RV data to a Keplerian orbital solution. The calculated un¬ 
certainties of the parameters account for both the RV and the 
stellar mass uncertainties. In this section we study five eclipsing 
binaries (Group A), six fast rotating KOIs (Group B), and pro¬ 
vide some hints for the presence of stellar-mass companions to 
two yet unsolved KOIs in our sample (Group C). We note that 
the nature of the transit on KOIs of both Group B and C remains 
unsolved with the current data. In the former, due to the impossi¬ 
bility of getting accurate RV measurements and the later because 
of the lack of good phase coverage. Tables 0] and [4] summarize 
our RV analysis for the three groups of KOIs. We also review 
the high-resolution ima ges of these KOIs already published in 
[Lillo-Box et alJ d2014al) . 

3.1. False positives: undiluted eclipsing binaries (Group A) 

Among the candidates being presented here, we have found five 
candidates that show large RV variations, incompatible with the 
presence of planetary-mass companions. These variations are 
modulated with the same period detected by the transit method 
using Kepler data. In the following, we analyze each of these 
false positives. In Table[3]we present the fitted orbital and phys¬ 
ical parameters for those KOIs with a clear Keplerian fit to the 

6 This pipeline was developed by Sebastian Sanchez. 

7 This algorithm has been already used to confirm the planetary na¬ 
ture of Kepler-9 lb injLillo^Boxetaljj^OlT^H) and as a sp ectral energy 
distribution fitter in Rivi ere-Marichala r et all j2Q 131 , 1 20 141) . 

s http://bass2000.obspm.fr/solar_spect.php 
9 http://exoplanets.org/code/ 


RV. FigureQjshows the measured radial velocities and the corre¬ 
sponding RV models of these KOIsQ 

KOI-0340.01 

This planetary candidate was ann ounc ed by the Kepler team 
in the first release of the mission ( IBorucki et al.ll2010 ). and still 
remains as a planetary candidate in the CFOR Its Kepler light 
curve presents transit-like dims every 23.673188 ± 0.000014 
days. Our analysis of the Kepler light-curve of this KOI in Sect. 
§|4]provides the detection of a secondary eclipse. Hence, we can 
keep the calculated period for the comnanioirl. 

This object was already included in the cata l og o f eclips¬ 
ing binarie s*^ presented an d updated in lMatiievic et al.l (120121) . 
Also, ISanterne et alJ (1201 2l) used the SOPHIE instrument to ob¬ 
tain two RV measurements at quadrature phases assuming circu¬ 
lar orbit, and derived a semi-amplitude of K — 34.577 ± 0.074 
km/s with such assumption, concluding that this system should 
be a single-line spectroscopic binary. Our CAFE data (obtained 
at eight different phases) provides a highly eccentric orbit (e = 
0.513 ± 0.005) for the companion (in good agreement with the 
location of the secondary eclipse, see Sect. §|4}. The fitted RV 
semi-amplitude provides a value of K — 16.28 ± 0.11 km/s. 
Assuming a host mass of M* = 1.10 ± 0.07 M 0 provided by 
the CFOP, we derived a stellar mass for the transiting object of 
Mi sin i = 0.214 ± 0.006 M 0 . The L-S periodogram indicates a 
period for the RV data in good agreement to the period obtained 
with the transits method thus indicating that both effects are pro¬ 
duced by the same object (see Fig.Q]). Hence, the non-planetary 
nature of the object transiting KOI-0340 is clear from this anal¬ 
ysis so that we can firmly and definitively discard this candidate 
as a planet. 

KOI-0371.01 

A periodic transit signal of P or b = 498.3928 ± 0.0042 days 
was detected by Kepler for this star. No significant difference 
in the depth of the two single eclipses is found in this case, so 
we maintain the period of this KOI as calculated by the Ke¬ 
pler team. The analysis of our RV data (with a time span of 
722 days) shows large relative variations for the stellar host. 
We find a clear (although broad) peak in the RV periodogram 
corresponding to the same periodicity as the transit signal (see 
Rg. m. According to our fit of the RV modulations, we find 
that the companion source should have a minimum mass of 
Mi sin i = 0.445 ± 0.014 M 0 . The physically bound stellar com¬ 
panion is calculated to be at a - 1.424 ± 0.047 AU, with a highly 
eccentric orbit of e - 0.407 ± 0.008. Thus, we can reject this 
planetary candidate and establish its binary nature. 

KOI-0686.01 

This planet candidate was released in the second catalog of Ke¬ 
pler planet candidates. The RV analysis reveals large variations 
with a clear modulation. The analysis of the AstraLux images 
of this object (see Sect. § EU >, yielded a 0.1% of probability for 
this KOI having a non-detected blended background star. Thus, 

10 KOI-3725 is not presented in this figure because its large v sin i pre¬ 
vents an accurate RV study. 

11 Note that if differences in the odd/even transit depths are found we 
should double the period since we would be seeing the two eclipses of 
an eclipsing binary. 

12 http://archive.stsci.edu/kepler/eclipsing_binaries.html 
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Table 2. Kepler Objects of Interest analyzed in this work. Observational data obtained with CAFE at Calar Alto Observatory. 


KOI 

f«Kep 

RA 

DEC 

Period^ 

MJ b> 

Rcomp W 

Nspec 

T 

span 


(mag.) 

(HH:MM:SS.SS) 

(DD:MM:SS.SS) 

(days) 

(M 0 ) 

(Re) 


(days) 

Group A 

- Eclipsing binaries 







340.01 

13.1 

19:50:39.52 

47:48:05.01 

23.673188 ±0.000014 

1.10 ±0.07 

0.231 ±0.021+ 

21 

321.1 

371.01 

12.2 

19:58:42.28 

40:51:23.37 

498.3928 ± 0.0042 

1.55 ±0.15 

0.77 ±0.39 

29 

722.0 

686.01 

13.6 

19:47:21.78 

43:38:49.64 

52.513565 ±0.000015 

0.84 ±0.12 

0.1248 ±0.0038+ 

19 

39.0 

3725.01 

10.1 

19:40:48.69 

38:31:10.38 

1.5704970 ±0.0000004 (rf) 

2.52 ± 0.53 

0.825 ± 0.007 

10 

17.0 

3853.01 

10.6 

19:11:13.73 

43:11:19.62 

21.512905 ±0.000050 

1.55 ±0.38 

0.31 ±0.22 

25 

70.9 

Group B 

- Fast rotators 







12.01 

11.4 

19:49:48.90 

41:00:39.56 

17.855149 ±0.000010 

1.31 ±0.27 

0.122 ±0.060 

47 

114.0 

131.01 

13.8 

19:56:23.41 

43:29:51.32 

5.014233 ± 0.000002 

1.31 ±0.26 

0.082 ± 0.029 

29 

418.0 

366.01 

11.7 

19:26:39.40 

38:37:09.32 

75.11202 ±0.00010 

1.69 ±0.17 

0.097 ± 0.042 

44 

72.9 

625.01 

13.6 

19:06:15.31 

39:32:04.09 

38.13819 ±0.00021 

1.46 ±0.24 

0.43 ± 0.21 

52 

66.0 

972.01 

9.3 

18:48:00.07 

48:32:32.00 

13.118908 ±0.000095 

2.02 ± 0.07 

0.055 ±0.013 

51 

335.3 

3728.01 

12.3 

19:35:05.31 

44:38:18.49 

5.5460830 ± 0.0000013 

2.05 ±0.17 

0.1880 ±0.0009 

40 

243.0 

Group C 

- Unsolved cases 







1463.01 

12.3 

19:13:02.07 

43:22:35.27 

1064.2681410 ±0.0001434 

1.17 ±0.28 

0.149 

14 

691.2 

3890.01 

13.2 

19:09:52.29 

37:57:59.90 

152.82630 ±0.00099 

1.76 ± 0.34 

0.3213 ± 0.0091 

22 

305.2 


Notes. (a> Transit periods are obtained from lBorucki et al1d201 ll) and jBatalha et alj d2013l ).Stellar masses are obtained from lHuber et al.l d2014l ). 
<c> Radius of the transiting compnions according to lBatalhaetalj j[^l_3|iJThis radius was obtained from our own transit fitting (see Sect. ffl.tln 
this case, the radius of the companion was obtained from iDiaz et alj d2014l) . (rf) In Sect. § [4] we have demonstrated that the actual period of this 
system is twice this value, i.e., 3.1409940 ± 0.0000008 days. 


it is highly probable that these RV variations come from a bond 
companion to the KOI-0686 system. 

The analysis of the detected transits by the Kepler team pro¬ 
vided an orbital period of P oy b = 52.513565 ± 0.000015 days. 
If we assume the same period and mid-transit time for the RV 
signal, we can obtain a model solution for the RV data. This fit 
provides a stellar mass for the companion transiting object of 
M b sin i - 0.089 ± 0.018 M 0 , thus confirming its stellar nature. 
We can thus confirm that KOI-068 6 is an eclipsingbinary. We 
note that a contemporary work by IPfaz et ah d2014l) also ob¬ 
tained the same conclusion, with the fitted parameters in good 
agreement within the uncertainties. 

KOI-3725.01 

The observations for this object were started while its status 
was "not-dispositioned" in the Kepler catalog. This meant that 
the light-curve analysis had not yet passed all the required cri¬ 
teria to be considered as a planet candidate. After the target 
was observed with CAFE, the Kepler team updated its status 
to "false positive". Our analysis of the light-curve in Sect. § [4] 
shows a clear difference in the depth of the odd/even eclipses, 
thus tagging this system as an eclipsing binary. Given this, we 
must consider a double period than the published by the Ke¬ 
pler team (i.e., P or b = 3.1409940 ± 0.0000008 days instead of 
1.5704970 ± 0.0000004 days). 

Our high-spectral resolution data (and light-curve analysis in 
Sect. [4]» confirms that this system is an eclipsing binary, showing 
large RV variations over the observed time span. However, the 
small number of phases acquired and the contaminated spectrum 
by the presence of two possible sets of lines, being one very 
broad ( suggesting large rotational velocities), prevents us from 
performing any detailed RV analysis. 

KOI-3853.01 

This KOI has been reported as a possible false positive in 
Article number, page 4 ofQ! 


the Kepler database^ The derived period according to the light 
curve analysis is / J or b = 21.512905 + 0.000050 days. The analy¬ 
sis of our high-resolution spectra obtained in 14 different phases 
reveals a large RV variation. When fitting the RV with a single 
object, an important periodic trend is found in the residuals. We 
tried to fit the RV data by assuming the presence of a third body 
in the system (see Fig. |3. The fit in this case was significantly 
improved. We can measure the significance of such improvement 
by using the Bayesian Inference CriterioiQ (BIC). We obtain 
a difference in the BIC value of ABIC > 100, showing a very 
strong evidence for the presence of a third companion. 

Interestingly, the period of the less massive companion (i.e., 
component C) is similar (in agreement within 3cr uncertainties) 
to the more massive companion (i.e., component B), thus being 
in a near 1:1 mean motion resonance (MMR). The periodogram 
of the RV data shows a double (blended) peak around the pe¬ 
riod corresponding to the the eclipse observed in the Kepler data. 
Both peaks have a false alarm probability (FAP) below 1%. In or¬ 
der to test the validity of the second periodicity in the data, we 
fixed the fitted parameters obtained for component B and we ran 
10 x 10 5 MCMC chains to fit the residuals of the data. We let 
the period of the possible C component to vary in the range 1-71 
days. The posterior distribution of the combined chains clusters 
around the estimated period provided by the joint fit (within the 
uncertainies). 

Thus all these three tests performed to our RV data (BIC, pe¬ 
riodogram, and MCMC analysis) favor the existence of a third 
component in 1:1 MMR against the two-bodies only configura¬ 
tion. The orbits of the two possible companions would be ori¬ 
ented in opposition, i.e., the difference in the argument of the 
periastron of both eccentric orbits agrees with being 180°. Ac¬ 
cording to the fitted models, and assuming the stellar mass for 

13 http://exoplanetarchive.ipac.caltech.edu 

14 A BIC = (x\ — X 2 ) ~ {ki ~ ki) hi N, where where ^ and k\ are the 
X~ statistic and the number of free parameters for the one-companion 
model and^ 2 and k 2 are the^y 2 statistic and the number of free parame¬ 
ters for the two-companions model. N is the number of RV points. 
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the A component provided in Table [2j the minimum masses of 
the two companions would be Mg sin is = 0.527 ± 0.044 M 0 and 
Me sin ic = 0.030 ± 0.016 M 0 . Note that the mass difference be¬ 
tween both components is significantly large, with the less mass 
component in the substellar regime. 

The unusual architecture of the proposed scenario would re¬ 
quire a more detailed study of its stability (which is beyond the 
scope of the present paper). Simulations of 1:1 me a n moti on 
resona nces in the planetary-mass regime by lAntoniadou et alJ 
( 2014 1) showed that these configurations are stable under spe¬ 
cific mass ratio conditions, mc/mg < 0.0205. Although our ob¬ 
jects are more massive, we note that this condition is fulfilled 
(within the uncertainties) for the KOI-3853 system. However, as 
stated before, more work is needed to analyze the stability of this 
system. 

Stellar pulsations induced by the companion could also ex¬ 
plain the correlated residuals in the RV fit. As we will see in 
Sect. [4] the light-curve of this target presents the so-called heart¬ 
beat effect, produced by tidal int eractio ns induced by the more 
massive stellar companion. In [Willems & Aertsl (120021) . the au¬ 
thors describe and calculate the RV variations induced by tides 
in close binary systems. According to them, there are specific 
resonant configurations that can lead to radial velocity variations 
even larger than 5 km/s, thus being a relevant explanation for the 
presence of the modulated residual of this close binary. How¬ 
ever, an accurate description of this system is out of the scope 
of this work. Future analysis studying the heartbeat-like modu¬ 
lation and the stellar modes induced by the close massive com¬ 
panion could unveil or reject the existence of the substellar com¬ 
panion in the 1:1 resonance orbit. 

3.2. Fast rotators (Group B). Upper limits to companion 
masses 

Planets orbiting fast rotating stars are more difficult to detect by 
the RV method. The line broadening produced by this effect pre¬ 
vents an accurate measurement of the line shift, increasing sig¬ 
nificantly the uncertainties for large v sin (i) values. Six fast ro¬ 
tators (v sin (i) > 20 km/s) were found in our sample. We com¬ 
puted individual spectral templates with the same physical prop¬ 
erties (effective temperature, surface gravity, metallicity, and ro¬ 
tational velocity) as the KOI by using the ATLASOT^l Then 
GAbox was used to compare the template and the observed spec¬ 
tra and determine the shift between both (i.e., the RV), by follow¬ 
ing the theoretical equivalency between maxim um-likelihood 
and cross-correlation presented in IZuckeJ (I2003H 1 6 1 

As stated, we can just provide an upper limit to the mass 
of the companion of these KOIs. We have assumed circular or¬ 
bits for all of them , as well as the lack of blended compan¬ 
ions contaminating the RV data (see Table|T}. A simple circular 
model was fitted to the RV data by fixing the period and time of 
mid-transit calculated by the Kepler team from the analysis of 
the eclipse. Thus, only two free parameters were fitted, the RV 
semi-amplitude ( K ) and the systemic velocity (V sys ). We ran one 
MCMC chain of 10 6 steps and discarded the first 10% to avoid 
dependency with the priors. The results for all six fast rotators 
are presented in Table [4] According to these results, all KOIs 
but KOI-3728 are compatible with having K - 0 km/s (non 
detection) within the uncertainties. Thus, we can take ctrv — 
K + 3 x ctk (where ctk is the uncertainty in K) as the upper limit 

15 http://kuruez.harvard.edu/grids.html 

16 See lLillo-Box et afl d2014bl) for a complete description of the method 
with GAbox. 


Table 4. Upper limits for the mass of the transiting companion to KOIs 
with large v sin i values (fast rotators, Group B). 


Group B ■ 

■ Fast rotators 




KOI 

v sin i 

K 

(km/s) 

Vsys 

(km/s) 

(TK 

(km/s) 

*Cp sini 
(M Jup ) 

12.01 

70 ±5° 

0.140 ± 1.425 

— 18.6 ± 1.1 

1.565 

25.2 ±3.7 

131.01 

26 t - 4 

1.400 ± 1.995 

-9.0 ± 1.8 

3.395 

35.8 ±5.0 

366.01 

35 ± 2 C 

0.027 ± 0.255 

8.5 ±0.2 

0.282 

8.70 ± 0.59 

625.01 

25 ± 2 C 

0.120 ±0.510 

-25.87 ± 0.42 

0.630 

14.1 ± 1.6 

972.01 

120 ± 5 C 

2.200 ±3.300 

-14.1 ±2.4 

5.500 

106.8 ± 2.6 

3728.01 

55t-“ 

4.810 ±2.775 

-52.6 ± 2.4 

7.585 

111.6 ±6.3 


Notes. * No uncertainties provided. <a) Values obtained from the analy¬ 
sis of TrES spectra by A. Byerla. lh> Values obtained from the analysis 
of Lick spectra by H. Isaacson/ Values obtained from the analysis of 
McDonald or TrES spectra by S. Quinn. 

of the semi-amplitude of the RV variations caused by a non- 
detected companion. In the case of KOI-3728, the solution for 
K is not compatible with zero, thus suggesting a possible detec¬ 
tion of the RV variations provoked by the transiting companion. 
According to this analysis, the companion would have a mass of 
0.068+0.041 M q . However, due to the large uncertainties and for 
consistency with the structure of this paper, we prefer to main¬ 
tain this KOI in this group B, although highlighting the possible 
detection of a stellar-like companion. The large v sin i, instead, 
does not allow us to ensure such detection. 

Hence, the RV amplitude ( K ) induced by the companion (if 
any) detected in these KOIs should be smaller than the afore¬ 
mentioned upper limit (cr S y), i.e. K < ctry- This amplitude is 
defined for circular orbits as 

2ttG( M c 3 omp sin 3 (0 \ 

P (Mcomp + M +) 2 J ’ U 

where the term in parenthesis is called the mass function ( Fm ), 
M* is the stellar mass, M comp is the mass of the companion, and 
P its orbital period. We can now determine the maximum mass 
of the companion that accomplishes K = (r KV . This value rep¬ 
resents the maximum mass (M™“ p sin i) of the hypothetic com¬ 
panion orbiting the star and inducing the detected transit signal. 
The calculated maximum masses of the hypothetic companions 
C^comp s * n 0 ^ listed in Table[4] 

All but two KOIs (KOI-0972 and KOI-3728) have upper 
mass limits below the substellar regime (~ 80 Mj up ). Indeed, in 
one case (namely KOI-0366), the upper mass limit is below the 
deuterium burning limit (~ 13 M j up ), suggesting that the object 
transiting this star is likely a planet. 

3.3. Unsolved cases (probable false positives, Group C) 

KOI-1463.01 

Only four phases (14 spectra in total) have been obtained with 
CAFE for this long-period candidate (_P or b = 1064.2681410 ± 
0.0001434 days). 

Our RV data (with a timespan of 691 days) shows large vari¬ 
ations of ~ 7.5 km/s, not corresponding to a planetary-mass ob¬ 
ject (see Fig. [3] left panel). However, the low number of epochs 
acquired and the timespan (shorter than the orbital period) pre¬ 
vents us from extracting definite conclusions about this system. 
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Table 3. Results of the fitting to the observed radial velocities. 


KOI 

K 

e 

CO 

\ sys 

M„ rap sin i<°> 

inc.< 4 > 

Fit® 

a 

F-test^ 


(km/s) 


(deg) 

(km/s) 

(Me) 

(deg.) 

(Mj up ) 

(AU) 


Group A ■ 

■ Eclipsing binaries 








0340 B 

16.28±0.11 

0.5138±0.0059 

239.71±0.58 

-84.007±0.045 

0.2144±0.0058 

89.95 

7.000±0.091 

0.1614±0.0054 

> 100 

0371 B 

8.252±0.086 

0.4070±0.0081 

94.75±0.47 

-56.855±0.081 

0.445±0.014 

89.95 

23.15±0.82 

1.424±0.047 

> 100 

0686 B 

6.3±1.7 

0.52±0.15 

63± 12 

-31.30±0.65 

0.089±0.018 

89.38 

0.88±0.45 

0.259±0.015 

> 100 

3853 B 

26.98±0.48 

0.3956±0.0084 

130.8±1.3 

-75.19±0.15 

0.527±0.043 

50.42 

35.5±2.2 

0.175±0.014 

> 100 

3853 C 

2.1±1.3 

0.56±0.20 

312±21 

-75.19±0.15 

0.030±0.016 

- 

0.011±0.019 

0.171±0.014 

> 100 

Group C ■ 

■ Unsolved cases 









1463 B 

> 3.3 

>0.18 

[80,200] 

[-28.5,-25.5] 

>0.183 ±0.030 

89.95 

>4.15 

2.15 ±0.19 

> 100 

3890 B 

> 2.5 

> 0.33 

[0,82] 

[-28.7,-24.2] 

> 0.097 ± 0.014 

89.95 

> 0.27 

0.676 ± 0.040 

3.6 


Notes. <a) Companion mass assuming primary mass from iHuber et~ahl (12014 1. see Table [2] <b) Orbital inclination provided by the Kepler team. 
(c) Mass function, defined as: F M = M/ omp sin 3 i/(M comp + M*) 2 . This parameter is directly obatined from the RV fit and does not depend on the 
primary mass. ^ F-test comparing the significance of the RV model to a straight line. 


We have used a radial velocity model including non-zero eccen¬ 
tricity to fit the available data. The free parameters used were 
thus the RV semi-amplitude, the eccentricity, the argument of 
the periastron, and the systemic velocity of the system. 

We ran a Monte-Carlo Markov-Chain with 10 6 steps to esti¬ 
mate these parameters and detect the correlations between them 
due to the lack of a complete sampling of the orbital phases. The 
degeneracy in the parameter space is important but we must re¬ 
mark that all solutions provide eccentricities larger than e > 0.18 
and semi-amplitude values larger than K > 3.3 km/s. This corre¬ 
sponds to a lower-mass limit for the companion transiting object 
of M”™ = 0.183 +0.030 M q . In the left panel of Fig. [3 we show 
a random sample of solutions obtained for this system. 

Thus, the transiting companion, according to the present data 
is likely of stellar nature, suggesting that KOI-1463.01 is a false 
positive. The solution to this system is still open and more RV 
data in a larger time span (since Kepler will no longer observe 
this target) will help solving this dichotomy, unveiling the na¬ 
ture and properties of the transiting object. We also note that 
iHirano et ali (120121) suggested that this candidate could actually 
be a late-type stellar companion transiting the star due to their 
analysis of the flux variations. 

KOI-3890.01 

This planet candidate was detected to transit its host star with a 
periodicity of 152.82630 + 0.00099 days. Our RV analysis, in a 
time span of 305.2 days shows a clear trend with a variation of 
3 km/s (see Fig. [3| right panel). Since we do not detect transit 
depth differences in the odd/even transits, we can fix the period 
and mid-transit time to that measured in the light-curve analysis 
by the Kepler team. Also, due to the small range of phases cov¬ 
ered by our observations (<f> e [0.6,0.8]), we cannot constrain 
the eccentricity and argument of the periastron of the orbit. We 
then proceed in the same manner as with KOI-1463, running an 
MCMC chain and analyzing the range of parameters. By doing 
so, we obtain a minimum RV semi-amplitude of K > 2.5 km/s. 
This translates into a minimum mass for the transiting compan¬ 
ion such that Mb > 0.097 ± 0.014 M 0 , with a minimum eccen¬ 
tricity ofe > 0.33. In the right panel of Fig. [3] we showarandom 
sample of solutions for this degenerated system. 

However, we must warn that in this case, the F-test compar¬ 
ing the RV model to a simple straight line provides a value of 
3.6, slightly favoring the false positive detection against a null 
detection model. Thus, although this candidate is likely a bi¬ 
nary, we cannot extract definite conclusions given the current 


data (mainly due to the small phase coverage). More RV at dif¬ 
ferent phases is necessary to better understand this system. 

4. Light-curve analysis 

In this section we analyze the Kepler light-curves of the targets 
studied in this work to look for possible interesting features that 
can contribute to better characterize these systems (secondary 
eclipses, photometric natural modulations - ellipsoidal, beam¬ 
ing, or reflexion-, and signatures of pos sible additional objects ) 
in the manner as we did in Kepler-91 dLillo-Box et al.ll2014d) . 
We note that the accuracy of some parameters (such as the ec¬ 
centricity) could be benefited from a joint analysis of the RV and 
the LC. However, in the particular cases studied in this paper, the 
eccentricity of the orbits are well constrained with the RV data 
(since most of the orbits are well sampled) and we assume the 
period and time of transit from the analysis of the LC by the Ke¬ 
pler team. Thus, we consider that the joint analysis would not 
improve the results significantly in the present cases. 

In the next subsections, we individually analyze the results 
of the targets were additional effects (rather than the detected 
transit by Kepler) have been found. 

KOI-0340 

We have detected the secondary eclipse (see Fig. ED of this 
RV-identified binary orbiting in a highly eccentric orbit (see 
Sect. § 13.Ik We have used a simple fitting model of a binary 
box with linear ingress and egress dependencies (enough for the 
purposes of the paper) to measure its depth, location, and du¬ 
ration. The period and time of mid-eclipse was assumed from 
that calculated by the Kepler team. The results show a depth 
of 62 - 650.6 + 6.6 ppm, located at <p — 0.31862 + 0.00002, 
and lasting d - 0.01005 + 0.00023 in phase units, equivalent to 
d - 5.70 + 0.13 hours assuming the period provided by Kepler 
team. By assuming the primary eclipse depth provided by the 
Kepler team (<5 [ = 22335 ppm, no uncertainty provided), we can 
derive a surface temperature ratio with 


where c)) is the depth of the deepest eclipse, 5 2 is the depth of 
the smaller eclipse, Ta/Tb is the ratio of effective temperatures 
of both stars. By using this equation, we obtain a temperature 
ratio of Ta/Tb = 2.421 + 0.006. 
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Also, as expected due to the large eccentricity and orientation 
of the orbit (obtained by the RV analysis), the location of the 
secondary eclipse is not centered at mid-time phase but instead at 
0 = 0.31862 ± 0.00002. The expected phase difference is known 
to be given by (I Wallenquistll 1950b 


A <p = 0.5 + e cos u) 


1 + (1/ sin/) 2 


(3) 


According to the parameters obtained with the RV analysis 
we would have expected the secondary eclipse to be located at 
4> = 0.335 ± 0.010, if we assume the inclination provided by the 
Kepler team. The difference between the expected time of mid¬ 
eclipse and the measured time is 9.3 + 5.7 hours. The source of 
this discrepancy is unknown but could be due to an underestima¬ 
tion of the uncertainties in the RV parameters or an inaccurate 
determination of the time of eclipse, translated into a small shift 
in the phase-folded light curve. 

Due to the large eccentricity found for this system in 
Sect. 13.11 we decided to perform_a dedicated modeling of the 
primary eclips40. We used the iMandel & Agol (120021) light- 
curve models and obtained the quadratic limb-darkening pa¬ 
rameters by interp olating the stella r parameters provided by 
iHuber et~afl (1201 4) to the values in IClaret & Bloemenl (1201 it) 
for the Kepler band. We fixed the eccentricity and argument of 
the periastron to the values found in the radial velocity analysis 
and leave the inclination (i), semi-major axis to primary radius 
(ci/Ra) and radius ratio ( Ra/Rb ) as free parameters. The results 
for this fit are: i — 89.6° + 0.2°, ci/Ra — 21.82 + 0.34, and 
Rb/Ra - 0.1452 + 0.0007. This model is shown in the upper 
panel of Fig. [4] The main difference with the Kepler team pa¬ 
rameters 0' K ep = 89.95°, (ci/Ra) Ke P = 14.31983 + 0.00722, and 
(Rb/Ra) Ke P = 0.142771^ “®^) is the semi-major to stellar ra¬ 
dius ratio. This is due to their assumption of circular orbit, while 
we have found a relatively large value for the eccentricity from 
both the RV and the secondary eclipse analysis. 


KOI-3725 

The most remarkable imprint in the light-curve of KOI-3725 
when phase-folding it with the period provided by the CFOP is 
the clear presence of two sets of eclipses with different depths. If 
we phase-fold the LC with twice the current period, we obtain a 
deeper eclipse at phase cf> — 0 than at (f> — 0.5. This is a clear sign 
of an eclipsing binary. According to this, the actual periodicity 
of this binary system is P or b = 3.1409940 + 0.0000008. 

The temperature ratio between both stars can be estimated 
from the depth ratio of both eclipses following Eq. [2] In this 
case, a simple fitting to the gaussian-like eclipses (enough for the 
purposes of this paper) provides Si = 2732 ± 30 ppm and 62 = 
1803 ± 35 ppm (see fitted functions over plotted in the bottom 
panel of Fig. [5j. This implies a temperature ratio of Ta/Tb — 
1.109 ±0.006. 


shape corresponds to a highly inclined orbit with the periastron 
oriented behind the plane of the sky. An accurate model fitting of 
this modulation is beyond the scope of this paper. Future works 
will focus on this object and fully characterize it. 

However, the heartbeat is not the only imprint in the light- 
curve of this system. A transit-like shape appears in the middle of 
the tidal modulation in the phase-folded light-curve. The Kepler 
team identified this transit as a potential candidate being of plan¬ 
etary nature. In the lower panel of Fig. [6] we show a zoom to this 
region once the "heartbeat" effect has been removed by a simple 
cubic spline fitting. The eclipse has a maximum depth of around 
2000 ppm. However, it is filled by some data points at lower 
depths. The period of both the heartbeat signal and the eclipse 
are exactly the same, indicating that both effects are provoked 
by the same object (i.e., companion B). Interestingly, we can see 
both a duration and depth variation in the individual eclipses (see 
Fig-0. Figure |8]displays the individual (detrended) eclipses. 

Our RV analysis (see Sect. 13.11 suggests that the motion 
variations of the primary star (A) are well described with the 
presence of two companions, a massive stellar object (compo¬ 
nent B) and a sub-stellar companion (component C) in a near 
1:1 mean motion resonance. Both objects would revolve in ec¬ 
centric orbits around the primary (massive) star. Based on this, 
we here suggest some possible explanations for the depth and 
duration variation of the detected eclipses: 

- Orbital wobble due to satellite companion. The presence of 
a third body orbiting component B can periodically vary its 
impact parameter, thus changing both the duration and the 
depth of the eclipse. Since we know from the shape of the 
heartbeat effect that the orbital inclination should be high, 
this should be a grazing eclipse of component B over com¬ 
ponent A. In the first periods of Fig. [8] the projected location 
of component B in the plane of the sky would be at impact 
parameters/? > I+Rb/Ra so that we see a small or no eclipse. 
By contrast, in the last periods observed, the impact parame¬ 
ter decreases so that the eclipse is deeper. A longer time span 
could have confirmed this hypothesis by detecting periodic 
variations of the eclipse depth. 

- Orbital wobble due to component C. In our RV analysis we 
found that the observed RV signal is much better modeled 
with two co-orbital stars with the same periodicity (in 1:1 
mean motion resonance) but in two eccentric orbits in op¬ 
position, i.e., with their pericenters differing by 180°. Com¬ 
ponent C, maybe in a higher inclined orbit (since we do not 
see any other dip in the LC), could be perturbing the orbit of 
component B, provoking an orbital wobble able to change its 
impact parameter and thus the eclipse depth. 

- Orbital circularization. If the system is being circularized, 
with component B orbiting closer to component A every pe¬ 
riod, the eclipse depth and duration will be increased with 
time. Since within the time span of Kepler data we do not 
see a new turn back to smaller depths, we cannot discard this 
possibility. 


KOI-3853 

The light-curve of KOI-3853 shows an interesting effect just 
discovered in a few number of tidally distorted binaries, the 
heartbeat pulsation. The shape of the light-curve modulation 
du e to this effect is well understood and was parametrized by 
iKumar et al.l d 19951) . In particular, in the case of KOI-3853, its 

17 The Kepler team assumes circular orbits in the transit fitting. 


With the current data, none of the previous conclusions can 
be discarded so that the problem remains open. However, our 
RV analysis suggests the presence of a third component in near 
1:1 mean motion resonance. Future observations and study will 
bring more light and try to unveil the nature of this system. In 
particular, a dedicated stability study is required to check the dy¬ 
namical plausibility of the co-orbiting substellar object in 1:1 
mean motion resonance. As we mentioned in Sect. § 13.11 the 
mass ratio between the two companions would be in agreement 
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to the allowed ratio for a triple system in 1:1 mea n motion res¬ 
onance as calculated by lAntoniadou et ail (120141) . although we 
warn that this study was specifically performed for planetary- 
mass systems. 


5. Summary and discussion 

We have analyzed the RV data obtained with the CAFE spec¬ 
trograph for 13 KOIs and the Kepler light-curves of some of 
them. In the majority of the eclipsing binaries found in this work 
(Group A), we found that the lighter co mpo nent is a ve r y low- 
mass stars (VLMSs). IZhou et all (|20 1 4l) and iDfaz et all (120141) 
published the last results on this kind of objects (see also refer¬ 
ences there in), adding more examples to the reduced number of 
well-characterized VLMSs (with detected eclipses and RV varia¬ 
tions). Added to this, our RV analysis suggests substellar masses 
for the transiting companions of the fast rotating stars (Group B). 

In Fig. |9j we show the mass-radius diagram for the com¬ 
panions found around these KOIs. In the case of KOI-0340 B, 
the radius provided by the Kepler team was well below the ex¬ 
pected value provided by any isochrone. We then decided to 
make a dedicated primary eclipse fit (see Sect. ED accounting 
for the large eccentricity measured in Sect. 13.11 By assuming 
the primary radius obtained by iHuber et~ak ( 2014 ). the new ra¬ 
dius of the companion, although larger, still lies below the ex¬ 
pected value. This could be attributed to an underestimation of 
the radius of the primary star ( Ra ). We can estimate R A from 
the a/R A ratio obtained by the primary eclipse fitting and the a 
value obtained in the RV analysis. The resulting updated radius 
is Ra — 1.589 ± 0.021 R 0 . If we now use this value and the radii 
ratio between both components obtained from the eclipse fitting 
( Rb/Ra X we obtain Ra = 0.231 + 0.021 R e . This new value per¬ 
fectly agrees with the expected value for the calculated mass of 
the companion. 

In the case of the fast rotator KOI-0972, the calculated upper 
mass limit for its transiting objects (assigning this mass limit to 
every candidate) also lie below the isochrones of VLMSs. The 
Kepler team detected two planet candidates in this system, with 
measured radii 7?oi = 6.1 ± 1.4 7?® and R^i — 1.73 ± 0.57 R s . 
There are thus two explanations for the discrepancy in the lo¬ 
cation of these objects in the mass-radius diagram: i) an under¬ 
estimation of the stellar radius, or ii) a real planetary nature for 
the candidates, with corresponding masses well below our upper 
limits. The former explanation is plausible since the host seems 
to be a sub giant star and the Darmouth models used to determine 
the radius could be biased for this kind of stars (see iHuber et al.l 
120141 ). The latter possibility is represented in Fig. [9]by the shad¬ 
owed regions. We can see that they would agree with having 
planetary-like densities, in agreement with their measured radii, 
if their masses are well below our detection limits. 

According to the mass-radius diagram and isochrones, we 
can also obtain an estimation of the radius of the suggested third 
component in KOI-3853 (i.e., KOI-3853 C). By assuming the 
calculated mass obtained by RV and the isochrones between 1- 
5 Gyr, the radius of this object should be in the range Rr ~ 
0.08 - 0.11 R q . 

In the following sections we summarize the results of each 
group of KOIs: 

5.1. Eclipsing binaries (False positives , Group A) 

We have spectroscopically identified five eclipsing binaries in 
the sample of Kepler planet candidates via RV analysis. Among 


them, KOI-0340 was already identified as a binary and KOI- 
3725 was initially classified as "not-dispositioned" but then up¬ 
dated to "false positive" after we had already obtained the spec¬ 
tra. We detected the secondary eclipse of KOI-0340 using Kepler 
data. By measuring its depth, we concluded that the temperature 
ratio of both stars is T^/Tf^ = 2.45 ± 0.03. Also, the location of 
the secondary eclipse agrees (except for a small offset) with the 
high eccentricity of the orbit derived by the RV analysis. In the 
case of KOI-3725, we concluded that the period of the system is 
actually twice that published by the Kepler team, since we detect 
clear differences in the odd/even eclipse depths. The location of 
the secondary (less deep) eclipse confirms a circular orbit for this 
system. 

Apart from the already confirmed false positive KOI-O 686 , 
we have reported other two eclipsing binary systems in this work 
(namely KOI-0371.01 and KOI3853.01). They present RV vari¬ 
ations with the same periodicity as the transit signal. The case of 
KOI-3853 is of particular interest. Its light-curve shows the clear 
signs of a recently discovered class of tidally distorted binaries, 
also known as heratbeat stars. We have found that the eclipse 
detected by Kepler changes progressively in depth and duration. 
The high inclination of the orbit of the system (known because 
of the shape of the heartbeat effect) agrees with this dim being a 
grazing eclipse of the two stars. Our RV data suggests the pres¬ 
ence of a third (substellar) component in the system revolving 
in an eccentric orbit around the primary star in a near 1:1 mean 
motion resonance with the more massive detected companion. 
Taking this into account, we hypothesize three possible expla¬ 
nations for the depth and duration variations: i) a third body 
perturbing the impact parameter on every transit, ii) an orbital 
wobble of component B due to the perturbing substellar com¬ 
panion, and iii) a possible circularization of the orbit with the 
stars getting closer and so decreasing the impact parameter. This 
dichotomy could be solved with a larger photometric time span. 
In the first and second scenarios we should see a modulation of 
the eclipse depth while in the second we should just see a rise 
in the depth value followed by a stabilization when the circu¬ 
larization is reached. However, with the current data, it is not 
possible to discard any of these configurations. The second sce¬ 
nario would agree with the observed RV, but it would require 
a detailed stability analysis. Other configurations could also be 
possible but more data are needed. 


5.2. Fast rotators: upper mass limits to Kepler planetary 
candidates (Group B) 


We have presented upper mass limits for six candidates with high 
rotational velocities. Because of this, our CAFE observations 
were unable to detect any RV variations below particular lim¬ 
its. These limits were used to set upper limits on the masses of 
the possible companions. We found that the companions (if any) 
to four of these KOIs (namely KOI-0012, KOI-0131, KOI-0366, 
KOI-0625) have substellar masses (i.e., M < 0.08 M 0 ), one of 
them (KOI-0366) having a maximum mass below the deuterium 
burning limit and thus being likely a true planet. On the con¬ 
trary, we found that KOI-3728 is likely a false positive, since our 
analysis in Sect. § 13.21 suggests the detection of a RV variation 
compatible with a stellar mass. However, the large rotational ve¬ 
locities prevents us from establishing a definite nature in both 
systems. 
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5.3. Unsolved cases (probable false positives, Group C) 

We have analyzed two cases whose nature remains unsolved af¬ 
ter our RV analysis. Although both of them show large RV vari¬ 
ations that strongly suggest that they are binaries, we could not 
definitively unveil their nature because of the lack of enough 
phase coverage and the low number of data acquired. In both 
cases we have provided limits for the different parameters in the 
RV fit. 

In KOI-3890, we obtained a minimum mass for the tran¬ 
siting companion of Mb > 0.097 M 0 , with the eccentricity of 
the orbit being e > 0.33. In the case of KOI-1463, we obtain 
Mb > 0.183 M q and e > 0.18. Thus, in both cases, the current 
data strongly suggests that the transiting companions are likely 
of stellar nature. 
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KOI-0686.01 


KOI-3853.01 





Fig. 1. Radial velocity analysis of the detected binary systems. In each set of images for every KOI we show the periodogram of the radial velocity 
data (upper panel), the fit tot he data according to Sect. l3.1l (middle panel), and the residuals of the fit (lower panel). 
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Fig. 2. Radial velocity analysis of the detected components of KOI-3853. Upper panel shows the RV solution for the whole system (black solid line) 
and the two contributions of the B component (dashed blue) and C component (dashed purple). Lower panels show the individual contributions of 
component B (left panel) and C (right panel). 



Fig. 3. Radial velocity analysis of the unsolved systems (Group C, see Sect. l3.3l l. We show one hundred random solutions accepted in our MCMC 
run (see Sect. l3.3l >. The residuals to the median solution are found in the lower panel 
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0.2 0.4 0.6 

Phase 



Fig. 4. Dettrended Kepler light curve of KOI-0340. In all panels grey 
small circles represent the Kepler data and the red filled circles is a bin¬ 
ning of the whole data with a bin size of 0.005 in phase. Upper panel 
shows a zoom to the primary transit and the fitted primary eclipse with 
a solid black line (see Sect. 0 . Middle panel shows the detected sec¬ 
ondary eclipse and the fit obtained in Sect.[4]as a solid black line. The 
lower panel shows the residuals of the fit, having a standard deviation 
of 173 ppm in the original dataset and 36 ppm in the binned light-curve. 


Fig. 5. Dettrended Kepler light curve of KOI-3725 and phase-folded by 
using the new calculated period (i.e., P mb = 3.1409940 ± 0.0000008 
days), twice that published by the Kepler team. The unequal depth of 
odd/even transits is a clear evidence of the binary nature of this object. 
Upper panel shows the entire phase of the system. Lower panel shows a 
zoom to the odd and even transits (in different colors and symbols). The 
primary transit has been shifted to 0 = 0.5 for comparison purposes. 
The transit model is shown as grey solid lines. 



0.996 


Fig. 6. Dettrended Kepler light curve of KOI-3853 and phase-folded 
using the period published by the Kepler team. The upper panel shows 
the entire phase to show the heartbeat effect. The red line shows the 
spline function fitted to remove the effect. Vertical dotted lines show the 
region zoomed Fig. [8] 
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Phase 

Fig. 8. Individual eclipses in the KOI-3853 system discussed in Sect. §[4] The color-code represents the transit number (equivalent to the julian 
date) since BJD = 2454966.6 days. The filled circles show the location (timing and depth) of the peak of the eclipse for each individual transit. We 
have omitted time spans in which Kepler data is not available at the phases of interest. The light-curve has been firstly detrended and then removed 
from the heartbeat effect as explained in Sect. § [4] 



Fig. 9. Mass-radius diagram of the companion objects with M < 0.3 Af 0 
found to some of the analyzed KOIs in groups A and B. The open square 
shows the location of the possible third body in the KOI-3853 (i.e., KOI- 
3853 C), whith an ad hoc assumed radius of 0.1 R Q (the dotted vertical 
error bar indicates our unknown about its actual size, see Sect. §0. 
We have used the i soch rones of 1 Gyr and 5 Gyr (from top to bottom) 
from Baraffe et alJ Il99 8h for solar metalicity (dotted dark blue line) and 
[Fe/H]=-0.5 (dashed li ght blu e line). In the case of masses below 0.1 
M e , we have used the Baraffe et al. 1 2003 1 isochrones (dotted-dashed 
gray lines) of 1, 5, and 10 Gyr (from top to bottom). The arrows in 
the Group B KOIs indicate that their masses are just upper limits. The 
dashed regions illustrates the range of masses allowed for the KOI-0972 
planet candidates (see Sect. § 0- For reference, we have plotted as 
light blue plus signs the known planets with provided masses and radius 
published in the exoplanet encyclopedia (www.exoplanet.eu). 
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Table A.l. KOI-0012: measured radial velocities. The S/N in Tables A.l to A. 13 is calculated as the median signal over the measured scatter of a 
continuum region close to 5500 A. 











Julian Date 

S/N 

RV 

Julian Date 

S/N 

RV 

Julian Date 

S/N 

RV 

(days)-2456000 


(km/s) 

(days)-2456000 


(km/s) 

(day s)-2456000 


(km/s) 

66.442058 

14.4 

-18.3!“ 

72.463525 

15.5 

-18.1™ 

87.634276 

15.5 

—18.5™ 

66.454881 

24.0 

-9.4“;> 

76.518732 

18.1 

-19.2134 

90.630156 

21.5 

-18.9™ 

66.467705 

19.0 

-9 6 +5 6 

-5.7 

76.533549 

34.2 

-H.5™ 

115.379299 

21.0 


69.439270 

15.0 

-18.4™ 

76.548365 

28.8 


115.391505 

29.9 

-18.0™ 

69.451609 

21.1 

—23.0;™ 

77.602751 

22.5 

-21. 4 ™ 

115.403712 

21.1 

-20.ru 

69.463949 

14.7 

-19.5™ 

77.615244 

30.5 

-19-0™ 

117.381785 

22.3 

-18.7™ 

70.434806 

24.8 

- 1 8- 1 ™ 

77.627738 

20.5 

-13.5™ 

117.393940 

27.0 

-14-3™ 

70.447044 

31.4 

20.62J | 

78.457637 

23.3 

—19-3™ 

117.406095 

15.0 

-20.9/“ 

70.459281 

19.1 

-21.8™ 

78.472721 

29.5 

-19.0™ 

129.381270 

17.0 

-21.9/4* 

71.431319 

23.0 

-18.9™ 

78.487805 

17.8 

-14.8™ 

129.393473 

27.8 

-2O.O/2* 

71.456236 

25.8 

-19.2/[° 

80.406010 

19.8 

-21.3™ 

129.405676 

21.9 

-18.7™ 

71.532211 

53.8 

-16.6™ 

80.418276 

31.2 

-19.8™ 

180.442531 

21.3 

-21.2™ 

71.608528 

30.3 

-16.6™ 

80.430543 

24.1 

- 17 . 7 ™ 

180.464569 

23.0 

-19-5™ 

71.632761 

29.8 

-19.4 ™ 

82.602553 

20.3 

-19-6 ™ 

180.464653 

41.1 

-19.6™ 

72.439229 

23.2 

-17.9™ 

82.614673 

26.2 

-20-21';* 

180.486860 

25.9 

-20.2/J4 

72.451377 

27.8 

-17.7™ 

82.626794 

16.6 

- 16 - 7 -.:v 





Table A.2. KOI-0131: measured radial velocities 


Julian Date 

S/N 

RV 

Julian Date 

S/N 

RV 

Julian Date 

S/N 

RV 

(days)-2456000 


(km/s) 

(days)-2456000 


(km/s) 

(days)-2456000 


(km/s) 

111.576836 

7.0 

-6 ? +1 - 7 
—1.4 

120.636824 

6.5 

-5 3 +31 

-3.4 

497.517894 

4.2 

-8-7!^ 

115.458790 

9.1 

-7-8/!;? 

418.567422 

6.5 

-12.1™ 

514.364588 

4.6 

-8.1™ 

115.493531 

7.4 

-i 0.7/2* 

418.599802 

7.6 

-“■8S1 

514.401807 

11.2 

-21.5/“ 

115.493595 

14.0 

-H-6/t;f 

418.632175 

12.2 

-10.5™ 

514.402164 

7.0 

-8.7/“ 

115.528462 

8.0 

- 8 . 1 ™ 

436.501911 

7.1 

159.32™ 2 

514.438668 

5.2 

-2 2+6.2 
J -- 5 -3.6 

119.544606 

4.5 

-8-4!?:? 

436.538388 

6.5 

0.0/4* 

529.484818 

8.7 

-12.1™ 

120.530916 

7.0 

-8.7™ 

436.573185 

5.2 

-10.0 ™ 

529.521132 

12.7 

-13- 4/40 

120.566218 

6.2 

-5.0/44 

497.439970 

5.2 

-3-533 

529.521166 

7.2 

-6.1™ 

120.583937 

13.4 

-7.9™ 

497.475725 

6.2 

-5-3 Hi 

529.557411 

4.5 

-12.9/*;4 

120.601789 

6.2 

-8-^ 

497.477863 

9.5 

-11.4™ 





Table A.3. KOI-0340: measured radial velocities 


Julian Date 

S/N 

RV 

Julian Date 

S/N 

RV 

Julian Date 

S/N 

RV 

(days)-2456000 


(km/s) 

(days)-2456000 


(km/s) 

(days)-2456000 


(km/s) 

111.395672 

13.8 

-103 499+ 0 045 

iuj.+^^_ 0045 

116.633232 

9.5 

-72.401™™ 

127.489427 

15.8 

—83.806™°* 

111.431211 

12.8 

-103 559 +0 095 

iuj.jj?_ 0095 

119.579846 

11.6 

-73.268/““? 

133.386589 

12.5 

—97.467/™“ 

113.605011 

12.5 

—91 99Q+ 0075 

-0.075 

119.612757 

19.6 

79 179+0.055 
/j.i /z,_o o54 

133.404004 

15.2 

—97 4‘S7+0.° 85 
y / / _ 0 074 

113.621494 

17.8 

-91 001 + ° 071 

y 1 • wwl _o.058 

119.614455 

12.4 

-73.262/°;“? 

133.421419 

8.1 

—Q7 S9+ 010 
-0.055 

113.637977 

12.5 

-90 909 +0 058 

vu.^w^_ 0082 

119.643968 

9.2 

-73381™°? 

137.604302 

9.5 

-86.483/°;“| 

116.598403 

7.8 

79 91 9+0.082 
/z,.z,u_ 0074 

127.454588 

15.5 

-83.850™? 4 

137.621876 

14.1 

-86.268™°° 

116.615818 

12.6 

79 999+0.074 
/ J -0-059 

127.472007 

21.9 

-83.789/°“° 

137.639450 

9.9 

- 86 - 171 ™?* 


Table A.4. KOI-0366: measured radial velocities 


Julian Date 

S/N 

RV 

Julian Date 

S/N 

RV 

Julian Date 

S/N 

RV 

(day s)-2456000 


(km/s) 

(days)-2456000 


(km/s) 

(day s)-2456000 


(km/s) 

66.494606 

15.4 

6.8/];? 

82.469641 

15.4 

8-3!?i 

112.644079 

15.4 

7 45+0.52 
' -0.46 

69.490469 

15.4 

"7 "7+1.8 
/- —1.6 

82.487003 

15.4 

7 Q4+0.78 
,y *-0.39 

118.374379 

15.4 

8 19 +0 - 47 

0,1 *-0.31 

71.558576 

15.4 

7 o <+0.37 
-0.28 

82.504365 

15.4 

7 77+0.84 
'•-0.25 

118.398555 

15.4 

0 70+O.4O 

0,/ -0.24 

71.570735 

15.4 

7 i c+0.38 
-0.31 

92.591566 

15.4 

q A7+0.40 
-0.53 

118.398689 

15.4 

9.00/SSS 

71.582894 

15.4 

7 50 +0 - 40 

/-JU -0.27 

92.615777 

15.4 

Q 1 1 +0.38 
y ‘ 1 —0.71 

118.423132 

15.4 

8.06/S;** 

72.599880 

15.4 

906™S 

92.615814 

15.4 

9 06 +0 - 33 

*- uu -0.30 

126.369758 

15.4 

"7 q-7+0.51 

7 0.61 

72.612031 

15.4 

8-80/S;*? 

92.640099 

15.4 

8.93™? 

126.382150 

15.4 

7 07+0.55 
-0.28 

72.624182 

15.4 

8.37/S;*S 

101.543893 

15.4 

7 17+0.44 
'-0.41 

126.394542 

15.4 

9 36 +0 - 64 
JU —0.26 

76.440887 

15.4 

9 00 + °- 93 

-0.42 

101.561229 

15.4 

8.64/S;*? 

130.380236 

15.4 

7 Q5+0.36 
' -0.51 

76.466954 

15.4 

Q <4+0.65 

101.578566 

15.4 

8-69/SJS 

130.392427 

15.4 

8 29 +0 - 31 
°- z, *-0.42 

76.467735 

15.4 

Q qa+0.88 
7-7 -0.29 

110.385278 

15.4 

7 07+0.62 
,yj -0.39 

130.404619 

15.4 

9 46+0- 32 

*•^-0.24 

76.492240 

15.4 

8.27/°°° 

110.397526 

15.4 

7 46+0-65 
' ' hu -0.56 

139.381895 

15.4 

7 7C+U.57 
' '' -0.51 

79.543631 

15.4 


110.409774 

15.4 

7 <7+0.88 
' -0.32 

139.394228 

15.4 

8-35/SI? 

79.555913 

15.4 

8-67/°™ 

112.619832 

15.4 

8.67/S;*? 

139.406562 

15.4 

8 -4i!S:S 

79.568194 

15.4 


112.631955 

15.4 

7 - 36 -047 
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Table A.5. KOI-0371: measured radial velocities 











Julian Date 

S/N 

RV 

Julian Date 

S/N 

RV 

Julian Date 

S/N 

RV 

(day s)-2456000 


(km/s) 

(days)-2456000 


(km/s) 

(days)-2456000 


(km/s) 

121.454974 

11.8 

-65.1691™“ 

505.601217 

6.5 

—53 219 +0049 

jj.2,17_o 053 

803.627383 

9.6 

-65 601 +O 055 

uj.uui_oo59 

121.467317 

17.7 

-65-13833S 

505.612199 

9.8 

cl i ci +0.049 
-0.059 

824.620384 

5.5 

-64 962 +0 067 
UH-.yuz.Q Q48 

121.479659 

13.2 

-65 200 + ° 053 
- 0 . 061 

505.623181 

7.5 

—57 q^+O - 058 

0059 

824.642337 

5.2 

-64 965 +0 070 

436.604913 

13.7 

-58 579+ 0 061 
uo.u / 7_ 0 _068 

528.419975 

10.2 

-53 651 +0 053 
jj.uji_oo49 

838.397849 

9.2 

-64 515 +0 066 
uh.jij_o.053 

436.629188 

13.7 

-58 399+ 0 068 
JO. 0 5 4 

528.438254 

18.5 

-53 526 +0054 
jj.jz-Uqoso 

838.416113 

14.1 

-64 386 +0 062 
UH-.JOU_0.049 

497.381256 

9.6 

-53 264 +0054 
jj.zuh_q 053 

528.456533 

15.7 

Cl C-7Q+0.059 
JJ.J/J_0.055 

838.434376 

10.2 

-64 406 +0 051 

UH.HUU_0 054 

497.394856 

15.2 

-53 194+ 0053 
JJ - i -0. 072 

559.495749 

11.7 

-58 079 +0 059 

JO.U/7_0 067 

843.450061 

11.3 

-64 246 +0 059 
UH-.Z.H-UQ048 

497.408455 

10.8 

—52 Q7 < s +0 - 053 

J 7 J -0.060 

560.496212 

14.5 

co oQ"2 +0.049 

38.ZOJ_o.070 

843.466246 

16.1 

-64 224 +0 053 
UH-.zz.H_o 062 

500.645274 

5.6 

m qco+U.072 

-0.049 

803.605301 

9.3 

-65.603!™™ 

843.482431 

11.2 

-64 270 +0 048 

UH.Z/U_0 053 

500.669675 

9.0 

-53 041 +0 060 

jj.uHi_ nn£i8 

803.616342 

13.5 

-65.618SSS? 





Table A.6. KOI-0625: measured radial velocities 


Julian Date 

S/N 

RV 

Julian Date 

S/N 

RV 

Julian Date 

S/N 

RV 

(days)-2456000 


(km/s) 

(days)-2456000 


(km/s) 

(days)-2456000 


(km/s) 

71.489198 

11.6 

-25 97 + 0 - 55 
"■*'-0.45 

88.573506 

1.7 

jo Q+4.2 
Z °'°-0.89 

94.589328 

9.8 

-25.77™°® 

71.506591 

16.5 

-25.80!°;“ 

88.592027 

9.8 

-26 27 +0 - 74 

zu.z / _q 67 

102.432371 

9.3 

-24 9 +11 

-0.56 

71.523984 

11.6 

-26 13 +0 - 65 
zu. i j_q 64 

88.597524 

13.0 

19 0 + 1.9 

J -0.99 

102.469438 

9.1 

-26 92 +0 " 
zu.7Z_o 54 

72.496141 

10.7 

-26.39!;™ 

88.627038 

7.8 

-24.6!!;’ 

102.469462 

15.6 

-31 2 +2 - 5 

J1 ^-0.72 

72.513475 

15.2 

_9< 9Q+0.84 
zJ'Z“_o.55 

91.459181 

8.3 

-25.9!^ 5 

102.506577 

7.8 

-24 3 +1 - 8 

^ -0.45 

72.530808 

10.7 

-25 5 +1 - 3 

J -0.49 

91.493991 

15.7 

-25 1 +L3 

" -0.75 

121.559629 

8.1 

-25 2 +u 

" -0.73 

77.443828 

11.6 

-25 57 +0 - 69 
ZJ.J/_Q8 9 

91.494063 

9.0 

-26 1 +1 - 2 
^°- a -0.84 

121.594577 

13.7 

-24 1 +L2 

zh,1 -3.8 

77.462844 

16.5 

-26 31 +0 - 59 
zu.ji_q62 

91.528730 

9.5 

-26 64 +0 - 95 
zu.uHq 90 

121.594628 

7.6 

-26.5™j 

77.481860 

11.8 

-26.66™;f 

92.525138 

11.2 

-25 37 +0 - 91 

zj.j »_Q 95 

121.629474 

7.8 

-24.7!!;® 

79.476176 

9.0 

-25 5 +LI 

J -0.92 

92.542705 

16.1 

-26.08™“ 

127.526771 

9.7 

-25.03!™° 

79.493735 

14.0 

-27.39!™ 6 

92.560272 

10.9 

-25.77™f 

127.561701 

9.0 

-26.1!!° 

79.511293 

10.4 

-25.08™“ 

93.453444 

10.2 

-25.0!};’ 

127.561732 

16.5 

-25.92™;®° 

80.462852 

8.2 

-29.7+1* 

93.471006 

14.6 

-24 2 +1 - 7 

zh.Z-i 2 

127.596724 

8.9 

-26 22 +0 - 61 
zu.zz_q 92 

80.487557 

13.2 


93.488568 

10.3 

-27 6 +L7 

137.461469 

8.5 

-25 l +u 

80.512261 

10.1 

-25.4™;° 

94.519592 

10.9 

-26 06 +0 - 94 
zu.uu _ 0 gs 

137.497308 

7.9 

-23.3![;° 

87.567872 

7.5 

-25 5 +4A 
ZJ-J- 1.2 

94.554529 

18.6 

-25 0 +L1 

—0.41 

137.498740 

15.2 

-27 1 +19 

z,/ -0.72 

87.585564 

10.0 

-25 7 +2 - 8 

^-0.95 

94.554666 

10.3 

-25.29™f 

137.537443 

9.5 

-26 52 +0 - 64 

ZU.JZ-1 oo 

87.603257 

6.2 

-24 7 +2 - 7 
_-0.71_ 








Table A.7. KOI-0686: measured radial velocities 


Julian Date 

S/N 

RV 

Julian Date 

S/N 

RV 

Julian Date 

S/N 

RV 

(day s)-2456000 


(km/s) 

(days)-2456000 


(km/s) 

(days)-2456000 


(km/s) 

91.564471 

9.5 

—29 i5't+ 0081 

zz7.i JJ_ 00 8j 

110.542846 

7.5 

-34 89 +012 

jh.o7_oo93 

124.524128 

6.5 

-34 625 +0 075 
jh.uzJq 069 

91.581828 

13.9 

-29 146 +0085 

z?.ihu_ 0 082 

110.560305 

11.0 

-34 77+ 011 

J -0.061 

129.435272 

10.8 

-33 343+ 0 - 081 

j j. jhj_ 0 085 

91.599185 

7.0 

_9Q '2Q7+0.082 
Z7. J7 /_ 0 073 

110.577764 

6.4 

-35 15 +on 
jj. i j_0 075 

130.540005 

10.0 

33 9'I7+0-0 7 6 
jj.zj ' _ 0 060 

100.433076 

8.2 

—24 7Q2 + 0 073 
zh./^z_ 0Q67 

124.453660 

6.9 

-34 620 +0 093 
jh.uzUq og! 

130.574986 

18.8 

-33.098™;“’ 

100.468077 

16.2 

—24 74Q+0.060 
zh. i*y _ 0 066 

124.488780 

8.2 

-34 500 + ° 061 
jh.juu_ 0076 

130.575012 

11.1 

-3 3 291 +0077 
jj.z^ioogg 

100.468144 

9.7 

-24.886™*’ 

124.488856 

13.1 

—34 S1 ^+0-068 
-0.077 

130.609940 

10.6 

-33 1SQ+ 0069 

100.503013 

9.5 

-24.834™* 6 
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Table A.8. KOI-0972: measured radial velocities 


Julian Date 

S/N 

RV 

Julian Date 

S/N 

RV 

Julian Date 

S/N 

RV 

(day s)-2456000 


(km/s) 

(days)-2456000 


(km/s) 

(days)-2456000 


(km/s) 

100.376405 

58.2 

-8-233 

114.503592 

53.2 

-8-533 

126.348892 

37.8 

-8.8!“ 

100.388716 

86.3 

-12-9!“ 

114.651427 

40.6 

-10.0!“ 

127.358229 

47.4 

-n.2!“ 

100.401026 

64.0 

-15-412 

115.359156 

38.5 

-n.3!“ 

128.359248 

34.8 

-14.9!“ 

101.610586 

41.2 

-12-812 

115.507564 

44.8 

-14.1!“ 

129.361016 

37.8 

-15 5 +3 - 3 

101.622675 

61.7 

-11 7 +4 - 5 

11 -5.8 

115.655973 

22.7 

-10.7!“ 

130.359776 

54.4 

-16. l“;= 

101.634765 

45.9 

-15.4!“ 

116.430815 

41.5 

-17.0!“ 

131.358377 

41.9 

_J4 1+3.8 

102.370604 

61.3 

-18.8!“ 

117.360460 

34.2 

-15.8!“ 

132.358859 

40.2 

_4+H 

—14 

102.382723 

89.3 

-9 3 +2 - 3 

y ‘ J -4.8 

117.399590 

41.8 

-13 1 +49 

AJ ‘ -5.3 

133.361742 

46.1 

-14 o +4 - 5 

A ^‘ -0.92 

102.394841 

65.1 

—1 7 +2 - 7 

'•'-4.5 

117.438719 

24.0 

_14 0+6.1 
A ^‘ —3.0 

137.356202 

41.2 

-22 5 +6 - 5 

103.602984 

62.3 

-13.2!“ 

118.353116 

34.0 

-18.9!“ 

138.359861 

43.1 

-19.0!“ 

103.615092 

87.4 

-16.412 

119.361354 

23.1 

-18-7 !=;= 

140.360137 

37.0 

7 7+5.0 
'•'-2.8 

103.627200 

61.2 

-10-712 

120.384560 

29.7 

-22.0!“ 

141.354705 

39.0 

-16.0!“ 

110.363523 

34.3 

-13.41= 

121.363682 

40.7 

-16.2“;= 

418.660632 

37.0 

-15.0!“ 

111.366741 

39.7 

-15.0!“ 

122.427018 

41.9 

-13.0!“ 

420.646028 

15.4 

77+20 

z '-5.9 

112.358116 

48.0 

-12-211 

123.355871 

30.5 

-6.7!“ 

420.658684 

21.2 

-82.22!“,“= 

113.357987 

34.7 

-6 4+ 5 - 9 
°‘ h -4.8 

124.352637 

41.6 

-14 2 +3 - 9 

-3.3 

426.552329 

13.7 

6+ 16 

-3.1 

114.355757 

34.3 

-11-812 

125.446545 

41.2 

-12 2 +3 - 5 
_ 1 -8.5 

435.646207 

34.0 

4+12 

-4.2 


Table A.9. KOI-1463: measured radial velocities 


Julian Date 

S/N 

RV 

Julian Date 

S/N 

RV 

Julian Date 

S/N 

RV 

(days)-2456000 


(km/s) 

(days)-2456000 


(km/s) 

(days)-2456000 


(km/s) 

114.381820 

12.8 

-30 76 +0 47 
ou. /u_ 047 

517.458750 

13.2 

31 76 +0 ’ 28 

J1 -' u -0.41 

838.477732 

12.2 

306.20!““ 

114.399160 

19.8 

-30 93+ 0 - 39 

036 

805.527896 

7.3 

305.37!““ 

838.494471 

16.1 

305.52!““ 

114.416499 

15.2 

-31.11!“? 

805.560558 

3.2 

306.09!“‘ 

838.511211 

8.9 

305.54!“;“ 

517.378382 

14.5 

31-18!g;i7 

805.560620 

11.3 

305.58!g;“ 

841.587718 

6.9 

305.43!“;““ 

517.418566 

20.1 

31.36!“ 9 

805.593404 

7.5 

305.57!““ 





Table A.10. KOI-3725: measured radial velocities 


Julian Date 

S/N 

RV 

Julian Date 

S/N 

RV 

Julian Date 

S/N 

RV 

(day s)-2456000 


(km/s) 

(days)-2456000 


(km/s) 

(days)-2456000 


(km/s) 

500.625109 

14.2 

-51.3!=;“ 

504.571903 

29.6 

-30 9 +1 - 2 

-0.97 

841.530484 

14.6 

-!5.0!“ 

501.383460 

19.2 

-46 2 +2 0 

-1.4 

504.609225 

22.8 

-41 9 +u 

* 1 - y -0.73 

841.536254 

20.6 

-29.°!“ 

503.425012 

23.2 

-42 1 +L4 

517.650133 

21.0 

-45.20!““ 

841.542025 

14.2 

-29.!!“ 

504.534581 

19.0 

-50 2 +L4 
_ -1.2 








Table A.ll. KOI-3728: measured radial velocities 


Julian Date 

S/N 

RV 

Julian Date 

S/N 

RV 

Julian Date 

S/N 

RV 

(days)-2456000 


(km/s) 

(days)-2456000 


(km/s) 

(day s)-2456000 


(km/s) 

354.369018 

10.1 

-36.3!=;= 

511.418066 

9.2 

_44+io 
^ -3.6 

523.540289 

15.3 

-48.6!=;= 

354.379993 

14.7 

-37.0!=;* 

512.472273 

6.9 

-54!“ 

525.518022 

14.4 

-47.9“= 

354.390968 

10.1 

-34 5 + ^- 3 

513.408866 

11.6 

-48.6!“ 

525.536150 

20.8 

—55.7“= 

498.570987 

12.2 

-57.8!“ 

513.427128 

18.5 

-50.7!“ 

525.554278 

14.8 

-53.9!“ 

498.583268 

16.9 

-58.4!“ 

513.445391 

13.7 

-50 9 +6 -° 

-6.9 

528.567258 

11.1 

-45 1 +4 - 2 

498.595550 

11.5 

-58.3!= 7 “ 

514.605639 

14.4 

-56 9 +5 2 

JU ‘*-4.3 

559.339189 

16.7 

-57.0!=;= 

500.576958 

4.8 

-30!=“ 

514.623762 

20.1 

-49““ 

559.371563 

15.8 

-57.1!=;“ 

500.602103 

6.8 

_67+ 25 

-3.6 

514.641885 

14.0 

-59.8!“ 

559.403933 

14.8 

-53 7 +5 - 5 

' -6.7 

501.428000 

3.9 

_49+ 17 

* y -8.2 

517.590081 

14.0 

-45 7 +4 - 9 

567.330058 

14.7 

-50 1 +26 

JU -4.8 

510.446736 

6.1 

_51+ 19 

J1 -9.7 

517.606266 

18.5 

-48.7!=;= 

567.367052 

15.0 

-47.4!=;= 

510.468900 

8.9 

-51+ 15 

J1 -4.4 

517.622450 

11.4 

-50.4“;“ 

567.403480 

14.6 

-50.2!=;= 

510.491065 

6.8 

_24+ 28 

^-5.7 

523.503712 

13.5 

-47 3 +3 - 7 

h '- j -5.6 

597.418633 

12.6 

-48!=| 

511.390550 

11.0 

-53 3 +4 - 5 

-4.9 

523.522001 

20.9 

-49.4“= 

597.418633 

9.3 

-47!^ 

511.404308 

14.7 

-56 3 +4 0 
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Table A.12. KOI-3853: measured radial velocities 


Julian Date 

S/N 

RV 

Julian Date 

S/N 

RV 

Julian Date 

S/N 

RV 

(days)-2456000 


(km/s) 

(days)-2456000 


(km/s) 

(days)-2456000 


(km/s) 

498.550774 

19.3 

-88.037!™? 

511.624458 

11.8 

— 57.982!»$4 

528.494410 

28.8 

_59 5 - 39 + 0.052 
-0.049 

500.427256 

6.8 

-78 182!°$: 

511.656257 

9.6 

-57.901!°™ 

528.515065 

43.6 

-59 420 +0 053 

503.522251 

23.8 

-67.622!°$7 

513.591688 

26.4 

-70 695 +0 057 

/u.u?j_ 0058 

528.535719 

32.7 

-59 298 +0 054 

J7.i7O_ 00 j3 

504.628523 

24.9 

-64 49 1+°°53 
uh-.h-?i_ 0053 

515.538671 

25.8 

-101 989 +0058 

iui.7O7_ 0053 

559.438923 

28.5 

-108.981!°™ 

505.425120 

25.4 

-62 606 +0 048 
uz.uuu _ 0 04g 

522.559858 

24.5 

-76 193+ 0053 
/U.17J_0 050 

559.449900 

39.6 

-108.981!°™ 

510.419882 

21.0 

-56 134+ 0052 
- 0. 052 

526.420908 

35.1 

-64 086 +0 050 
W,UOU - 0. 048 

559.460877 

27.4 

-109.009!°™ 

511.363906 

25.3 

C' 7 1 0^7 +0.051 
D/.15/_oo5i 

526.439128 

45.2 

-64 033 +0 048 

u+.ujj _ 0 0 5 2 

569.430575 

25.6 

-64 360 +0 053 

uh-. juu_ 00 53 

511.553563 

33.1 

S7 1 T "2 +0.056 
J/.1JJ_ 00 56 

526.457348 

27.8 

-63 921 +0 048 

UJ ^ zl - 0 . 0 5 4 

569.467568 

16.5 

-64 290 +0 057 

U+.Z?U_ 0048 

511.569631 

14.5 

-57 611 056 

_ J 11 -0.054 








Table A. 13. KOI-3890: measured radial velocities 











Julian Date 

S/N 

RV 

Julian Date 

S/N 

RV 

Julian Date 

S/N 

RV 

(days)-2456000 


(km/s) 

(days)-2456000 


(km/s) 

(days)-2456000 


(km/s) 

499.397905 

6.7 

— 27-368!°™ 

511.496364 

2.8 

-26 35 +016 

U-JJ -0.074 

516.552987 

9.4 

-25.955!°™ 

499.438485 

4.1 

-27.74!°g 

511.535899 

5.5 

-26 176 +0074 

z.u.1 /u_ 0052 

516.554295 

16.9 

R "3 c+0.064 
z.j.ojj _ 0 Q9 8 

502.629346 

6.2 

-27 330 +0 072 

513.484972 

7.6 

-26 418 +0 052 

z.u.h-io _ 0 072 

516.589283 

10.0 

-25.870!°™ 

510.526583 

5.0 

-26 60 +01 ° 
z,u.uv_ 0 100 

513.522807 

14.2 

-26 240 +0 055 

z,u.zh-u _ 0056 

522.447180 

3.2 

-24 7l +0 - 38 

1 -0.090 

510.564787 

2.5 

-26 35 +0 - 24 
z - u . JJ _0 2 4 

513.523404 

8.2 

-26 140 +0072 

^,0.1^ 0 065 

522.485030 

12.2 

-24 798 +0 096 
/ “°_ 0 055 

510.577098 

8.2 

-26.391!°™ 

513.560045 

8.5 

-26 260 +0 056 
z,u.zuu_ 007 3 

522.485601 

7.2 

-24 970 +0 098 

- 0.064 

510.639922 

7.1 

-26.387!°™ 

516.520616 

9.8 

—25 94?+0.065 
z,j.^z,_ 0066 

522.522308 

8.8 

-25 041 +0 070 

z.j.u+i _ 0096 

511.457622 

8.2 

-26.370!°?f 
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